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We study the spatiotemporal dynamics, in one and two spatial dimensions, of two complex fields which are 
the two components of a vector field satisfying a vector form of the complex Ginzburg-Landau equation. We 
find synchronization and generalized synchronization of the spatiotemporally chaotic dynamics. The two kinds of 
synchronization can coexist simultaneously in different regions of the space, and they are mediated by localized 
structures. A quantitative characterization of the degree of synchronization is given in terms mutual information 



measures. 
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Introduction spatial coupling, so that chaos is absent from them in 



I 1. 

O. 

^ Starting from the pioneering experimental study by Huy- 

^ gens with two marine pendulum clocks hanging from a 

L| common support |Huygens, 1665], synchronization phe- 

. ,—{ fiomena have been subject of intense study in many 

/\^ physical and biological systems |Winfree, 198C; Btrogatz 



H |and Steward, 1993| . Since [ Pecoraand Carroll, 1990| | es- 



spatially homogeneous situations. Spatiotemporal chaos 
has thus a rather different origin which could lead even- 
tually to a different kind of synchronization. The Com- 
plex Ginzburg-Landau (CGL) equation is one of such 
model systems: in the absence of spatial coupling, it is 
simply the normal form for a Hopf bifurcation, so that 
chaotic behavior is absent. But spatial coupling induces 



tablished that chaotic oscillators can also become syn- a huge variety of intricate chaotic behavior [shraiman 



[ihronized, many applications and extensions of the orig- 



inal idea have been identified. Some of them arc the pos- 
sibilitics of partial (i.e. phase) synchronization [[Roscn- 



blum et ai, 1996 1, generalized synchronization [Rulkov 



et ai, 1995; Kocarev and Parlitz, 1996a |, and synchro- 



et al 1992| ; phate, 1994a| ; phate, 1994b| ; [Eguiluz et ai, 
1998 f^ The possibility of synchronized chaos between 



nization of spatiotemporally chaotic systems | Amengual 



et ai, 1997|. A natural class of systems in which syn 



a pair of amplitudes satisfying a pair of coupled CGL 
equations in one spati al dimension was considered in 
[Amengual et ai, 1997 1. A kind of generalized synchro- 



chronization of spatiotemporal chaos can be explored is 



the one constr ucted in the following way 1 Kocarev and 
Parlitz, 1996b| : Take a couple of chaotic systems that 
synchronize when appropriately coupled. Then make 
copies of this composite system, one copy at each point 
of a spatial lattice, and couple them spatially to study 
how the spatial coupling modifies the synchronization 
characteristics. This class of systems displays a kind 
of spatiotemporal chaos, and of chaos synchronization, 
which is a natural extension of the one obtained from 
the single composite system. 

But there are spatially extended systems that dis- 
play a kind of spatiotemporal chaos produced by the 



nization was found and characterized. In this Paper we 
show that different kinds of synchronization are possi- 
ble, all mediated by the presence of localized objects in 
the equations solutions, which become specially robust 
in twodimensional situations. For small couplings, both 
usual and generalized synchronization coexist simulta- 
neously in different regions of the space. For larger cou- 
plings only generalized synchronization remains. 



2. Vector Complex Ginzburg-Landa^ = 0.7. in addition to the blue holes there also red 

maxima in the amplitudes. It is clear that the maxima 
in one of the components appear where the other com- 
ponent has a hole, so that a strong anticorrelation is 
present between these variables. 



equation in the one dimensional 
case 



In the context of amplitude equations, the CGL is the 
generic model describing slow modulations in the oscil- 
lations of s patially coupled osc illators close to a Hopf 
bifurcation | van Saarloos, 1994 1 . Pairs of coupled CGLs 
have been derived in a variety of contexts, mainly related 



to the interaction of counterpropagating waves | Amen- 
gual et iL, 1996 |. A different context in which coupled 
CGLs appear is in optics: the interaction of the two 
polarization states of light in large aperture lasers has 
been show to be described by coupled CGLs with par- 
ticular symmetries which allow the pair to be thought as 
a Vector Complex Ginzburg-Landau (VCGL) equation 
[3an Miguel, 1996]. The two components of the VCGL 
equation can be written as 

dtA+ = A+ + {1 + ia)\/^A+ 

-{l + ^(3){\A+\' + ^\A^\')A+ 
= A^ + {l + ia)V^A^ 

-il + ^P){\A^\'+J\A+\')A^. (1) 



dtA^ 



A+ and A 



are the two components of the vector com- 
plex field, which in optical applications are identified 
with the left and right circularly polarized components 
of light, a and /3 are real parameters. We will restrict 
our discussion to the case in which 7 < 1 is a real num- 
ber and 1 + a/3 > 1 is satisfied (Benjamin-Feir stable 
range). These restrictions appear for laser systems pre- 
ferring linearly polarized emission pan Miguel, 1996 1. 
Note that, different from other contexts in which cou- 
pled CGL Eqs. have been derived, in the optical po- 
larization context a group-velocity term is absent from 
(P. An extensive analysis of coupled CGL equations for 
other parameter regimes can be found in |van Hecke ei 
a/., 199|. 

Within the Benjamin-Feir stable range, there is a re- 
gion of parameters for which a single CGL Eq. has still 
a regime of spatiotemporally chaotic behavior. It is the 
so called spatiotemporal intermittency regime |Chate, 
1994a[|~^tates in this regime consist in patches of travel- 
ling waves interrupted by localized objects (depressions 
or holes in the field amplitude) that move rather errat- 
ically around the system while emitti ng waves and per- 
turbations. These homoclinic holes [ van Hecke, 1998 
are the responsible of sustained spatiotemporal chaos in 
the system. Both nonlinear dispersion (the (3 parame- 
ter) and spatial coupling are needed to obtain this kind 
of chaotic behaviour. Figure [^ shows a spatiotemporal 
configuration of the modulus of the field, obtained from 
one of the two equations in (|l|) when 7 = 0, so that 
it is uncoupled to the other component. Blue lines are 
the trajectories of the localized depressions disorganiz- 
ing the system, whereas green-yellow regions are lami- 
nar states. When two coupled equations with 7 9^ 0, are 
considered, new objects come into play. Fig. || shows 
the modulus of the two components |yl_|_| and |^_| for 




Figure 1: Spatiotemporal evolution of the modulus 
squared |A(a;,t)|^ of the scalar CGL equation (7 = 0) 
for a = 0.2 and /? = — 2. Time is running upwards from 
t — to t — 400 and x is in the horizontal direction. 
The system size is L = 512. 




Figure 2: Spatiotemporal evolution of \A^{x,t)\'^ (left) 
and |^_(a;,t)p (right) for 7 = 0.7. The other parame- 
ters are as in Fig. [^ 



When analyzed separately, both |yl+| and |yl_| dis- 
play spatiotemporal chaos. The knowledge of one of 
these variables, however, gives a large amount of infor- 
mation on the other, since they are strongly anticorre- 
lated. This is precisely the content of the concept of gen- 
cralized synchronization |Kocarev and Parlitz, 1996a; 
the two chaotically evolving vari- 



Rulkov et ai, 199! 



ables are not identical as in the usual synchronization 
case, but there is a functional relationship between them 
which allows close prediction of one of them when the 



other is known. In | Amengual et ai, 1997 1 the functional 



relation was identified and a mutual information calcu- 
lation showed that the generalized synchronization be- 
came more perfect as the parameter 7 approached 7=1 
from below. For 7 = 1 the functional relation between 
the two fields is \A+\'^ + |A_p = 1. The appearance of 
correlations between the components is clearly mediated 
by the presence of the localized objects, so that this is 
a kind of generalized synchronization of chaos specific 
of spatiotemporal chaos, that is, it is absent in systems 
without spatial dependence. 

It should be noticed that only the moduli |^+| and 
\A-\ develop correlations and we thus find amplitude 
synchronization. The corresponding phases do not be- 
come synchronized. This is the case opposite to the one 



commonly observed of phase synchronization [Rosen- 
blum et\ ai, 1996|. The reason for this is that the cou- 



pling between the fields in Eq. n involve only the moduli. 




mf% mm 
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Figure 3: Temporal evolution of |^+(a;,i)p (left) and 
|j4_(a;,i)p (right) starting from an initial condition 
where A+{x,0) « A_(a;,0). a = 0.6, /3 = -1.4 and 
7 = 0.7. The other parameters are as in Fig. H 



Although the meaning and origin of the generalized 
synchronization found is clear, one often prefers to apply 
the term "synchronization of chaos" to situations closer 



to the original ideas |Pecora and Carroll, 1990 1, which in 



the present context would mean a tendency of the two 
fields involved to take identical values, not just anticor- 
related ones. We have explored, numerically and analyt- 
ically, this possibility and found that there exist states 
in which the amplitude of the two moduli are identi- 
cal. They are however dynamically unstable, so that the 
system does not approach them unless particular initial 
conditions are carefully selected. One example is shown 
in Fig. |3[ during the early part of the evolution the two 
components are perfectly synchronized, but a decay to 
the previous anticorrelated case occurs at later times. 
Further states will be discussed elsewhere. Two types 
of localized structures can be seen in Fig. ||, the ones 
for which both fields have simultaneously a hole and the 
ones in which a hole in one of the fields is associated to a 
pulse in the other field. Hole-hole and hole-pulse local- 



ized structures will also be present in twodimensional 
systems as we will show in the next section. Further- 
more, we will see that topological restrictions can make 
stable in higher dimensions the hole-hole localized ob- 
jects related to the ones appearing at the early times in 
Fig. |. 



3. Twodimensional defects 



The onedimensional localized structures of the previous 
section become topological def ects in tw o dimensions. 
They have been introduced in |Gil, 1993 and properly 
classified in [ Pismen, 1992 ; Pismen, 1994[ . There are two 
main types of topological defects in (|lj) : vectorial defects 
are objects for which the two amplitudes become iden- 
tical near the defect core, where both vanish. They are 
the two dimensional analogs of the hole-hole structures 
present in the early part of Fig. ^. In the other class 
of defects only one of the amplitudes vanishes, whereas 
the other presents a maximum. These are the twodimen- 
sional analogues of the hole-pulse localized structures in 
Fig. gand in Fig. ||for larger times. Fig. |4|(a) shows the 
amplitude of one of the two fields containing both types 
of defects. The vectorial ones emit waves that entrain 
a whole domain around them whereas defects in which 
only one amplitudes vanishes behave more passively and 
remain at the domain borders. In fig. H(b) the global 
phase (f)g = 04. -|- (j)_ , where 0+ and (f>^ are the phases of 
A^ and A_, is plotted. Clearly, there are two different 
kinds of vectorial defects. One is produced by two de- 
fects of the same topological charge, which correspond 
to the two-armed spiral in the plot of 4'g. The other is 
generated with defects of opposite charges, leading to a 
target pattern in the plot of (j)g . 



In the dynamical evolution from random initial con- 
ditions, vectorial defects are spontaneously formed for 
the set of parameters o f Fig. H. This is the sa me set of 
parameters explored in [ Amengual et a/., 1997 1 for a ID 
problem and it leads here to a glassy or frozen configura- 
tion like the one shown. However, increasing the value 
of the coupling parameter 7, the vectorial defects be- 
come unstable. The system evolves then in a disordered 
dynamics which is governed by defects in which only 
one amplitudes vanishes. The number of defects is con- 
served, after a transient regime, during very long times 
of evolution. A snapshot of this state is shown in Fig |5| 
where it is seen that a zero of one of the amplitudes cor- 
responds to a maximum of the other one. These defects, 
strongly anticorrelated in amplitude, move together in 
time so that the kind of spatiotemporal chaos they sus- 
tain displays generalized synchronization between the 



two amplitudes \A^ 




For 7 ~ 0.3, the cloud of points broadens and the line 
\A+\ — \A-\ becomes diffuse and disappears. This be- 
havior identifies the instability of the vectorial defects. 
Such instability appears in a narrow range of 7 and with 
different mechanisms for the two types of vectorial de- 
fects. As the coupling is further increased (7 = 0.5 and 
0.95), the cloud of points approaches the curve given by 
jA+p -I- l^-P = 1 as in one dimension. This indicates 
anticorrelation between A^ and v4_ and a generalized 
amplitude synchronization of the dynamics which in- 
creases with 7. Similar qualitative behavior with 7 is 
observed for other values of a and /? for which a glassy 
state occurs at 7 = 0. 



Figure 4: Field configuration showing different kinds 
of defects in two spacial dimensions. One of the ampli- 
tudes, |A+p (a), and the global phase (j)g (b), at a given 
time are shown. Parameters used: 7 = 0.1, a = 0.2 and 
/3 = 2, system size 128x128. 




Figure 5: Snapshot of the intensities of the two fields 
(a) |A+|^ and (b) |A_|^ for 7 = 0.8 and a and P as in 
Fig 0. A dark blue (yellow) dot corresponds to a zero 
(maximum) of the amplitude. 







A quantitative description of the process of increas- 
ing generalized synchronization as 7 is increased can be 
given in terms of information measures as already intro- 



duced for the ID case in [Amengual et ai, 1997|. This 



description also identifies the transition from glassy to 
dynamic states due to the instability of the vectorial 
defects. The joint probability density p(|yl+|, |yl„|) is 
plotted in fig. || for different values of 7, as a 3D plot 
and, in the same figure, as a density plot. The density 
plot is obtained taking the simultaneous values of |^+| 
and \A^\ at different space-time points. For a small cou- 
pling, 7 = 0.1, we obtain a diffuse cloud of points with 
a broad maximum around \A^\ ~ |^-| — 1 with devi- 
ations from these values being uncorrelated, except for 
the points laying in the line |^+| = |^-|. The presence 
of such line in plots of coupled variables is the classical 
signature of conventional synchronization. The points 
on the line correspond to the core of the vectorial de- 
fects in which the two amplitudes take the same value. 



Figure 6: Joint probability distribution p(|A+[, [A_|), 
for a and /3 as in Fig Q, shown as a 3D surface for 
different values of 7. From top to bottom and from 
left to right, 7 = 0.1, 0.3, 0.5 and 0.95. On top of 
each surface the corresponding plots of |yl+(a;,i)| vs 
\A-{x,t)\ are shown (the density of dots is propor- 
tional to p(|A+|, |A_|)). A time average was taken over 
100 samples separated At — 1. The 3D surfaces for 
p(|A+|, |A_|) has been cut at p = 0.0025. 



Two quantities that can be extracted from the prob- 
ability density are the entropy for a single amplitude 
and their mutual information. The entropy is defined 
as H{X) = —^^p{x)h\p{x), where p{x) is the prob- 



ability that X takes the value x. H{X) measures the 
randomness of a discrete random variable X. For two 
random discrete variables X and Y, with joint proba- 
bility density p{x, y), the mutual information I{X, Y) = 
— J2x yP{^^ y) ^'^[p{^)p{y) I P{^ 1 y)] gives a measure of the 
statistical dependence between both variables; the mu- 
tual information being if and only if X and Y are 
independent. Considering the discretized values of |^+| 
and \A_\ at space-time points as random variables X — 
|A_|_|, Y — I A- 1, their mutual information is a measure 
of their generalized synchronization. The dependence 
of the entropy of ^+1 and \A^\ and their mutual infor- 
mation / on the coupling parameter 7, shown in Fig. 
M identifies three regimes: For low values of 7 there is 
a glassy or frozen configuration dominated by large is- 
lands around vectorial defects. This relatively ordered 
state gives a relatively small value of H . The mutual 
information is not large because of the weak coupling 
between the fields. A second regime is associated with 
the instability of the vectorial defects. This is high- 
lighted by a maximum value of the entropies and a min- 
imum of / for 7 ~ 0.3. This instability first disorders 
the configuration by reducing the size of the domains. 
This yields an increase of the entropies, but this disor- 
der is not correlated in both components as indicated 
by the decrease of /. For higher values of 7 (7 > 0.35) 
we enter the third regime characterized by disordered 
configurations which evolve in time like the one shown 
in Fig. ^ The disorder of each amplitude |^±| is mea- 
sured by a relatively large value of H , while the increas- 
ingly synchronized dynamics is measured by an increas- 
ing value of / which approaches its maximum possible 
value [/ = H{\A+\) = i?(|A_|)] as 7 ^ 1". 
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Figure 7: Entropy of 1^4+ 1 (squares) and |A_| (triangles) 
and their mutual information / (diamonds) as functions 
of the coupling parameter 7 for a and /? as in Fig 0. 



4. Conclusions 

We have described generalized synchronization phenom- 
ena in the spatiotemporal dynamics of two complex fields 



which are independent components of a vector com- 
plex field. Synchronized dynamics of the amplitudes 
occurs via localized objects. The motion of these struc- 
tures produce the dynamical disorder. In rf = 1 local- 
ized structures are always dynamical and have the form 
of anticorrelated pulse-hole structures. Hole-hole type 
structures are seen to decay after a short transient. In 
d = 2 the hole-hole type structures become topologi- 
cal vectorial defects which are stable for small coupling 
and produce frozen configurations. For coupling above a 
threshold value these vectorial defects disappear and the 
persistent dynamics is governed by structures which are 
topological defects of only one of the fields. These d — 2 
structures show strong amplitude-amplitude anticorre- 
lation and they are the analog of the d — 1 pulse-hole 
structures. A quantitative characterization of the de- 
gree of synchronization is given by mutual information 
measures. 
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